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ABSTRACT 


'An  analysis  of  generally  laminated  rectangular  com- 
posite plates  subject  to  lateral  pressure  and  hygrothermal 


loads  is  presented.  Included^xn  the  analysis1* are  the  ef- 
fects of  transverse  shear  and  transverse  normal  deformations—-  — 
✓tffie  latter  effect  being  included  to  determine  the  importance 
of  through  the  thickness  hygroscopic  expansions  of  the  poly- 
meric matrix.  Solutions  are  obtained  using  the  theorem  of 
minimum  potential  energy  for  both  clamped  and  simply-sup- 
ported boundary  conditions.  In  addition  various  in-plane 
boundary  conditions  are  studied  to  demonstrate  the  behavior 
of  unsymmetric  laminates  subjected  to  a hygrothermal  environ- 
ment. 


v7  Numerical  results  of  plate  deflections  and  rotations 
for  various  width  to  thickness  and  length  to  width  aspect 
ratios  show  that,  while  transverse  shear  deformations  are 
important,  transverse  normal  deformations  have  no  apprecia- 
ble influence  on  the  solution.  Results  presented  display 
the  magnitude  and  distribution  of  lamina  stresses  under 
lateral  pressure  and  hygrothermal  loads  for  different  in- 
plane restraint  conditions  and  laminate  lay-up  patterns. 


vi 


NOMENCLATURE 


parameter  in  moisture  diffusion  solution 


^ plate  length  in  x-direction 

A area  of  plate 

A:;  extensional  stiffnesses  of  laminated  plate 

8 plate  length  in  y-direction 


0. . bending-extensional  coupling  stiffnesses  of 
AJ  laminated  plate 


integration  constant 
D hygroscopic  diffusivity 


bending  stiffnesses  of  laminated  plate 


elastic  modulus  of  lamina 


tracing  constant  for  TNS 


body  forces 


stiffness  parameter  from  TNS 


shear  modulus  of  lamina 


vii 


viii 


lamina  thickness 


vectorial  distance  from  the  mid-plane,  z 
to  the  upper  surface  of  the  Kth  lamina 


H total  laminate  thickness 


h4  stiffness  parameter  from  TNS 
•rj  th  tt  -rHi 

a>  i,V i , a l parameters  defined  in  eqn.  (2.32,  2.33,  2.34) 

KJ.lC.K7 1 

VY\ parameter  in  moisture  diffusion  solution 

moisture  distribution  within  plate 

ka*  thermal  and  hygroscopic  moment  resultants, 

' 'a  i ' V respectively 

/T.nrr  parameters  defined  in  equation  (2.27) 

moment  resultant  normal  to  plate  edge 

moisture  concentration  at  upper  or  lower 
° plate  surface 

/VN.  moisture  concentration  at  upper  and  lower 
u'  L plate  surface,  respectively 

/ 

N number  of  lamina 


in-plane  stress  resultant  at  plate  edge 

T u 

N N:  in-plane  thermal  and  hygroscopic  stress  re- 

Aj  * sultant,  respectively 


parameter  in  diffusion  solution 


V 

*iYii 


yW  yW 

^ 3 


ft 


7«,? 

7w  if" 
© 

X, 

MS 

<5,% 

<*> 


strain  energy  density  function 

plate  coordinate  directions 

in-plane  rotations  in  the  x and  y direc- 
tions, respectively 

coefficients  of  thermal  expansion 

coefficients  of  hygroscopic  expansion 

transformed  coefficients  of  the  thermal  and 
hygroscopic  expansions  in  the  plate  coordi- 
nates, respectively 

convergence  test  parameter 

strain  components 

mid-surface  extensional  strains 

transverse  normal  deformation 

constants  based  upon  the  boundary  conditions 

angle  from  the  x*  to  the  1*  coordinate 
axis  for  each  lamina 

plate  curvature  strains 

Poisson's  ratio 

stress  components 


characteristic  beam  functions 


CHAPTER  ONE 


INTRODUCTION 


Investigations  into  the  use  of  reinforced  composite 
materials  with  polymeric  matrices  for  structural  applications 
have  documented  the  deleterious  effects  of  hygrothermal  en- 
vironments on  structural  performance  [Ref.  15,  16,  17,  18, ' 

19,  20,  21,  22].  The  deleterious  effects  consists  of  degra- 
dation of  the  matrix  dominated  mechanical  properties  and  hy- 
groscopic expansions  in  the  matrix  [16],  As  indicated  during 
the  AFOSR  workshop  at  the  University  of  Delaware  [16] , a need 
exists  for  a comprehensive  analysis  method  including  the  ef- 
fects of  a hygrothermal  environment. 

In  an  effort  to  partially  satisfy  this  need,  an  analy- 
sis of  generally  laminated  rectangular  composite  plates  sub- 
ject to  lateral  pressure  and  hygrothermal  loads  is  presented. 
Included  in  the  analysis  are  the  effects  of  transverse  shear 
and  transverse  normal  deformations.  The  latter  effect  being 
included  to  determine  the  importance  of  through  the  thickness 
hygroscopic  expansions  of  the  polymeric  matrix.  Solutions 
are  obtained  using  the  theorem  of  minimum  potential  energy 
for  both  clamped  and  simply-supported  boundary  conditions. 
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In  addition  various  in-plane  boundary  conditions  are  studied 
to  demonstrate  the  behavior  of  symmetric  and  unsymroetric 
laminates  subjected  to  a hygrothermal  environment. 

The  analysis  of  laminated  composite  material  plates 
with  various  boundary  conditions  has  been  developed  by  many 
researchers,  [2,  2,  5,  11,  12,  13,  23,  24,  25,  26,  33,  34, 
35].  Of  these  [3,  5,  26]  have  shown  the  importance  of  trans- 
verse shear  deformation,  particularly  for  highly  anisotropic 
materials.  Also  [24,  25]  include  thermal  expansion  loads  in 
their  analysis.  Only  [25]  addresses  the  problem  of  moisture 
expansion,  but  considers  moisture  and  thermal  expansions  sep- 
arately and  omits  the  effects  of  transverse  shear  and  trans- 
verse normal  deformations* 

In  regard  to  transverse  normal  strains  (TNS) , refer- 
ences [27,  28,  29]  have  included  TNS  to  investigate  the  norm- 
al stresses  generated  near  a free  laminate  edge.  Flaggs  [6] 
is  the  only  researcher  including  TNS  with  moisture  swelling 
in  an  investigation  of  the  stability  of  generally  laminated 
plates. 

The  analysis  of  hygroscopic  residual  stresses  gen- 
erated in  rectangular  plate  elements  without  boundary  condi- 
tions was  presented  by  Pipes,  Vinson  and  Chou  [1].  The 
treatment  of  the  hygroscopic  expansions  and  moisture  diffu- 
sion in  [1]  is  used  in  this  work.  However,  the  diffusion 
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<5* * 


where  C(j  » lamina  stresses  in  the  j direction, 

■ lamina  strains, 

yT* 

6:  » coefficient  of  thermal  expansion, 

v»* 

Jj^  * coefficient  of  hygroscopic  expansion, 

* temperature  distribution  through  the 
thickness  of  the  plate  at  time,  t, 
defined  as  the  difference  between 
the  temperature  at  that  point  at  time, 
t,  and  the  stress  free  temperature  of 
the  plate, 

= moisture  concentration  distribution 
through  the  thickness  of  the  plate 
at  time,  t,  defined  as  the  difference 
between  the  moisture  concentration 
at  that  point  at  time,  t,  and  the 
stress  free  moisture  concentration. 

Using  the  following  contracted  notation. 


* 

A-  X.tfi 

6aa  * 

<*] 

O'  _ 

6xi  * 6* 

oil  * 

6,i  - 6f 

^x  * 

*,x  * 

(2.2) 


---  ~ ..  r- 


ill  i 


f 


equation  (2.1)  can  be  expanded,  i.e 


The  Q_  are  the  elastic  stiffness  matrix  components  in  the 
material  coordinate  system.  These  are  given  in  Ref.  [2] 


for  the  three  dimensional  case  as 


Using  these  transformations,  the  generalized  stress-strain 


relation  for  the  Kth  lamina  becomes 


T(x,*>  - ftUzf) 
€y  - F/Tb,*)  - V W 
^-S'Ta.ti-Fx^rnM 


0.  Q*  G„  o o K)» 

5a  5a  5x3  o 0 

5*  o o zQ» 

o o O XQm^s  o 

0 o O 20^23%  o 
Q+  Q*  <5*  O o 
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The  transformed  stiffness  matrix  components  are  defined  as, 

(■ 

Qii*  Qh  m*  "V  2L  V y\  { Qhl  4 2-Q<*)  + h**  Oxx 
Qtx=  mV  (Qw  + Qml  “ H’Qtfc)  + (m*  trl* ) Qlx 

0.3*  mi*  0,3  4 a*  0a3 

Qfc  5 * m n3GU  + m3nQ«  ~ mn(w?-n'X^v+  £ Qt*) 

Qaa  * * On  ♦ 2. mV  (<?,x  ♦ cLQ*)  + **  Ozx 
0*3  s h*0is  + m*  Qx3 

033=  Q»  (2*10) 

= " m^nQxa.  + m a Qn  + mv>(  ma-  a1  )(Q,a  4 ZQu!) 

Q&  * “ "^($*3  “ Oa) 

Qhh  * m*  Osh  + AX  Qsc 
Oss-  s ~mn  (<?*h  " Qss) 

Osf  ~ vy'*Q«r  * ^ OfM 

Qfcb  - n (Q„  4 Qxx”  2»Q*X)  + (ma-  a1)  Q**  • 


The  transformed  coefficients  of  expansion  are, 

^ yT  x VT  ^ v^«  y«  \ 1 

• x ” ®ii  m + Tjx  yx  9: k **  »n  m •+  ^ a 

?yT  = r,7  n1  + j£  V"  vs1  ♦ 

s 2. { y,^ - Vj^wiv,  s *-h«- 


(2.11) 


Note  that  while  the  hygrothermal  expansion  coeffi- 
cients are  dilatational  in  the  material  coordinate  system, 
when  transformed  into  the  plate  coordinate  system,  coeffi- 


cients involving  the  in-plane  shear  strains  are  introduced. 


Strain-Displacement  Relations 

The  linear  strain-displacement  relations  in  the 
cartesian  coordinate  system  are. 


J Ur 
iy 

*»  ■ * if) 

(2.12) 

auj 

C S.  - 

€l  *» 

Only  the  linear  order  terms  are  included  in 
equation  (2.12)  in  that  only  small  deflections  relative  to 
the  plate  thickness  are  being  considered. 

The  assumed  form  of  the  u,  v,  w displacements  in  the 
x»  y,  z directions,  respectively,  are 

ir(x,Y,x}  - (2,13) 

UjCx.y.t)  s uj0(x,y)  ♦ { (z) 

where  the  subscript  "a"  denotes  the  mid-plane  displacements. 
The,  f,  in  the  w(x,y,a)  equation  is  a tracing  constant  for 
the  transverse  normal  deformation  function,  n(z).  Substi- 
tution of  equation  (2.13)  into  equation  (2.12)  gives. 


I 
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In  this  analysis  the  variation  of  u and  v with 

o o 

respect  to  the  z direction  will  be  assumed  zero,  therefore 


Application  of  the  Theorem  of  Minimum  Potential  Ener< 


The  solution  of  genera1  .ated  composite  plates 

under  uniform  pressure  and  hygrothermal  loads  will  be  obtained 
using  the  Theorem  of  Minimum  Potential  Energy. 


strain  energy  density  function 


where  W 


volume  of  the  elastic  body, 

ith  component  of  the  surface  traction 

ith  component  of  the  deformation, 


portion  of  the  surface  over  which  the 


The  mathematical  statement  of  the  minimum  potential  energy 


where  ^ ( ) represents  the  variation  operation. 

In  the  present  work,  body  forces  will  be  neglected.  The 


The  surface  tractions  considered  in  the  analysis  will  be 
those  normal  to  the  surface  of  the  plate,  specified  as 


Then 
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j T;  UA  <^S*  = JT  Pu,v)  dA*  (2.21) 

sJf  "a* 

In  addition,  the  plate  theory  assumptions  on  the  normal 
stresses  and  strains  will  not  be  used.  The  transverse  shear 
stresses  will  be  included  because  of  their  importance  in  the 
analysis  of  highly  anisotropic  plates. 

The  strain  energy  density  function,  W,  including 
the  hygrothermal  strains,  will  be 

W * 2[j°*  I1!  + €y'^T-  S'"* 


2*»[vKr- 


(2.22) 


Thus,  the  potential  energy  of  a generally  laminated  plate 
will  be, 

V = i & IjJ  [*[v  ?T‘  ^+<Jy[s-  *y  \ 


+ °i[«2 ' ?T  - fnj  + ^[2  ^ + 0-^2  e>J 
"■^[2^-  ^T-  ?"y  rnjj  dz.  d A* 


(2.23) 


Substitution  of  the  stress-strain  relations,  equation  (2.9), 
and  the  strain-displacement  relations,  equation  (2.16),  into 
equation  (2.23),  and  integrating  across  the  thickness  for  all 
terms  except  those  involving  n (z) , gives, 
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At  this  point  one  could  obtain  an  approximate 


solution  by  using  the  Rayleigh-Ritz  method.  This  is 
accomplished  by  assuming  a functional  form  of  the  unknown 
displacement  and  rotation  variables,  substituting  these 
into  the  potential  energy  equation,  and  taking  variations 
with  respect  to  the  unknown  amplitudes.  However,  problems 
arise,  in  that,  while  forms  for  the  uQ,  vQ,  wQ,  a,  and  6 
variables  can  be  assumed  for  certain  boundary  conditions, 
the  form  of  the  transverse  normal  deformation  function,  n» 
is  not  intuitively  known.  In  order  to  circumvent  this 
problem,  the  Euler-Lagrange  equation  for  the  6n  term  will 
be  used  to  approximate  n in  terms  of  the  other  displacements 
and  rotations.  Thus,  taking  the  variation  of  equation 
(2.24),  and  collecting  all  terms  involving  6n  gives,  after 
integration  by  parts. 


A *k-i 

[Sf'S?’ (£».'*  MpT 

*IC-I 

[ij;M5:?;*o^4t]nhji 


From  the  first  half  of  equation  (2.28)  either  0it  * o or. 


I 

I 

I 

[ 

B 

D 

B 

E 

B 

B 

|:B 

B 

E 

I 


20 

The  second  half  of  equation  (2.28)  represents  the  natural 
boundary  conditions  at  the  h = and  h^  edges  of  the 

plate. 

Equation  (2.29)  can  now  be  used  to  obtain  an 

approximate  expression  for  the  3n/9_  terms  in  equation  (2.24). 

z 

2 2 

Solving  equation  (2.29)  for  3 n/3_  and  integrating  with 

z 

respect  to  z gives. 


(2.30) 

+ ^[QaC  + •*  * §*  %y]T 

OW 

The  constant  of  integration,  C(x,y),  will  be 

arbitrarily  set  equal  to  zero,  making  z = 0 the  reference 

. 

temperature  or  moisture  concentration  for  the  effects  of 

' 

the  transverse  normal  strains.  Substituting  equation  (2.30) 
into  equation  (2.24)  gives, 


♦%[fc  jf 

+ Bit  17  4 tt ♦ *Efc  t + &$ 

+fUt  *8# + ff  $ 

♦AJ**  ♦-£**fck  +tlf]+A»tf  ■ ~t HHtfi 
<k[£  * e%  * tftyf]  - Mtiftf * Ift 

+ 8Jt  t + W + 1 1 + + fcftfljfr 
* * Mfcif  * $?fc] 


in- 


3 

g 

0 

1 

° 

0 


D 

1.2.6  jj 

1.2. 3. 6 

(2.33)  ij 

Q 

0 

0 

1.2.3. 6 |j 

(2.34) 

D 

0 

0 
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Equation  (2.31)  can  now  be  used  to  obtain  the 
solution  for  a generally  laminated  plate  by  the  Rayleigh- 
Ritz  method.  The  Rayleigh-Ritz  method  requires  that 
admissible  functions  be  chosen  for  the  unknown  displacements 
and  rotations.  These  functions,  which  satisfy  the  geometric 
boundary  conditions,  can  be  substituted  into  the  potential 
energy  equation  and  variations  taken  with  respect  to  the 
unknown  amplitudes.  This  leaves  a set  of  simultaneous 
equations  to  solve  for  the  unknown  amplitudes. 


* 


To  study  the  effects  of  the  transverse  normal  strain 
(TNS)  a tracing  constant,  f,  has  been  included  in  the 
derivation.  To  include  TNS  effects  set  f = 1.  To  exclude 


TNS  effects  set  f = 0. 


Hygrothermal  Effects 


The  effects  of  a hygrothermal  environment  on  a 
laminated  composite  plate  are  introduced  into  the  analysis 
through  dilatational  strains  in  the  constitutive  relations. 
These  strains  are  related  to  the  thermal  and  moisture  con- 
centration distributions,  T(z,t)  and  M(z,t),  by  the  coeffi- 
cients of  expansion,  as  expressed  in  equation  (2.3).  After 
integration  through  the  thickness  of  the  plate,  the  hygro- 
thermal strains  take  the  form  of  stress  resultants  and  couples 
as  defined  in  equation  (2.26).  To  evaluate  the  stress  re- 
sultants and  couples,  it  is  necessary  to  define  distributions 
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which  are  integrable. 

Transient  Case:  Pipes,  Vinson  and  Chou  [1]  derived  the 
moisture  concentration  distribution  using  the  classical 
diffusion  equation, 


d£M  _ An  = o. 

u*  a*  u 


(2.35) 


where  D is  the  diffusion  constant  of  the  material.  For 
symmetric  moisture  absorption,  the  solution  to  equation 
(2.35),  with  the  appropriate  boundary  conditions,  is. 


= /A0[  I - ^0'T'n  cos 


where, 

and 


(2.36) 
' (2.37) 


m4 


[hOL  eai°M 

0! 

U«*ie  ] 

t (2.38) 

1 

For  symmetric  moisture  desorption, 

— mC 

tAu.k)  « Mo  £ 

For  single  surface  moisture  absorption. 

Mm  = Mo[i-|ePnc0sV(x^)] 

where  - _ 


(2.39) 


(2.40) 
, (2.41) 


0 

D 

D 

0 

0 

0 
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p »■+)  tlL 

jfj  th+1 


(2.42) 


These  solutions  are  also  valid  for  the  transient 
thermal  distribution,  by  using  the  thermal  diffusivity  for 


Steady-State  Case:  The  steady-state  distributions  for  the 
temperature  and  moisture  concentration  will  be  simply  linear 
functions  in  the  z direction. 


/A(z)  - 


(2.43) 


(2.44) 


where  the  subscripts  u and  L signify  the  upper  and  lower 
surface  respectively. 


The  given  distribution  functions  M and  T can  be 
substituted  into  equation  (2.26)  and  integrated  to  obtain  the 
hygrothermal  stress  resultants  and  couples . 
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comings,  particularly  for  obtaining  the  stress  solution  for 
a problem.  The  convergence  of  the  deflection  solution  is 
generally  rapid,  assuming  that  a reasonable  displacement 
function  is  chosen.  However,  to  obtain  the  same  degree  of 
convergence  for  derivatives  of  the  displacement  functions 
that  are  needed  to  obtain  stresses,  more  terms  of  the 
admissible  functions  are  needed. 

Boundary  Conditions 

The  solutions  obtained  herein  were  for  plates  with 
the  following  three  support  configurations. 

1.  All  Edges  Simply  Supported 

2.  All  Edges  Clamped. 

3.  Two  Opposite  Edges  Simply  Supported, 

Two  Opposite  Edges  Clamped. 

In  any  analysis  of  unsymmetrically  laminated  plates, 
the  in-plane  displacements,  u and  v,  must  be  included. 

These  displacements  are  necessary  to  account  for  the  bending- 
extensional  coupling  present  in  unsymmetric  laminates. 

The  boundary  conditions  used  in  the  analysis  of 
generally  laminated  simply  supported  plates,  as  presented  by 
Whitney  111]/  are 

SI:  Wn  = 0,  Mn  = 0,  un  = 0,  Nnt  = 0 (3.1) 

S2:  Wn  = 0,  Mn  = 0,  Nn  ■ 0,  ufc  * 0,  (3.2) 
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where  the  subscripts  n and  t represent  the  normal  and 
tangential  direction  to  the  plate  edge,  respectively. 

For  the  clamped  plate  the  following  boundary  con- 
ditions are  used. 

Cl:  Wn  = 0,  an  = 0,  ufi  = 0,  Nnt  = 0 (3.3) 

C2:  Wn  = 0,  an  = 0,  Nn  = 0,  ut  = 0 (3.4) 

For  the  simple-clamped  plate  the  edges  parallel 
to  the  x axis  use  the  clamped  boundary  conditions.  The 
edges  parallel  to  the  y axis  use  the  simply -supported 
boundary  conditions. 

Assumed  Displacement  and  Rotation  Functions 

Contrary  to  the  approach  of  many  researchers 
[6,7,8],  which  include  transverse  shear  deformations,  the 
reduced  unit  width  beam  approximation  from  the  Euler- 
Lagrange  equations  will  not  be  used  to  remove  the  rotational 
dependent  variables,  a and  6.  Instead,  a functional  form 
for  these  variables  will  be  assumed  as  done  by  Wu  and  Vinson 
[9]. 

The  assumed  functional  forms  for  the  simply- supported 
boundary  condition,  SI,  are, 


. 
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GJCx.y)-  Z $ Sin  njy 

A 8 

«.(x,y)  = 1.  $ Q«  sw  *«<»»»¥ 

n%«V*i^  *»'AS  A ° 

/S(*,y)  « Jz  S.  CoSfi&*  SlV\nKX 

^ •'  rh*vi,f  **W  'Av'  A Is 

UM  = ^ £ "Uvhk  Sln^SJL  Cos  to* 

n-^r  K.i,»^  A B 

^,y)  = ^ ^ V^  Cc>s  ^ s»nJ5|y 


(3.5) 


(3.6) 


The  form  for  the  S2  boundary  condition  will  be  similar  to 
SI,  except  that  the  following  forms  for  u and  v are  used 

instead  of  equation  (3.6). 

oo  o*> 

U U,y)  = £ Z.  Uv..  coi  eaui  siU  toi* 

A 8 


(3.7) 


^Ay)'  £ V_  SintoSi  cos,  **y 

A O 


The  assumed  form  of  the  displacements  and  rotations 
for  the  Cl  clamped  boundary  condition  are  the  characteristic 
beam  functions,  originally  presented  by  Warburton  in  (41 
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and  extended  by  Wu  and  Vinson  in  [9], 


<U*>t,n<y) 


of>  QO 


(3.8) 


0 (x,y)  = 1.  £ A„k  <}Wy) 

U (Y,y)  r £ £ si^jeai*.  cos  jap 


©O  o-O 


v<* ,y)  = £ £ V, 


C©S  )2i32$  Slv\i 


(3.9) 


where 

W = cos  I>n  [j.  - -L1)  * cos V'/vfe-i) 

^nW  ' ^ i)  (3a<>) 

= Cos  + V 
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and 

% = -s jM/VzJ 

(/V./2-) 

+av\  (■/£»)  + WvV\  (/£aj)  a 0- 


(3.11) 


(3.12) 


By  interchanging  a and  6 and  substituting  y for  x,  n for  m, 
and  B for  A,  the  equations  given  in  (3.10)  yield  the  proper 
functions  for  the  y direction.  The  solution  to  the  transcen- 
dental equation  (3.12)  gives  the  necessary  values  for  the 
Um  constant  in  equations  (3.10).  These  solutions  are, 


yU,*  I .SOSbZ-K  , 2.HVV75X 

/“Ks  (*■■*•£)*■  Kz3 


Again  for  the  C2  boundary  conditions,  u and  v will  be 
changed  to  the  functions  given  in  equation  (3.7). 

It  should  be  noted  that  only  the  odd  mode  shapes  of 
m and  n are  involved  in  the  analysis.  This  occurs  because 
only  the  odd  terms  will  be  non- zero  for  the  loads  and 
symmetric  boundary  conditions  considered. 

For  the  simple-clamped  boundary  conditions  the 
following  combinations  of  the  previous  functional  forms 


are  used 
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The  u and  v variables  have  the  same  form  as  given  by  equa- 
tions (3.6)  or  (3.7)  depending  on  the  desired  in-plane 
boundary  conditions. 

Deflection  and  Stress  Calculations 

With  the  values  of  the  unknown  amplitude  coeffi- 
cients, W , T , A , U , and  V , the  deflections  and 

mn  mn  mn  mn  mn 

rotations  of  the  plate  can  be  easily  calculated  from  the 
appropriate  support  condition  given  in  equation  (3.5)  to 
(3.14). 

The  laminate  strains  can  be  calculated  using  the 
strain-displacement  relationships  in  equation  (2.16)  along 
with  the  assumed  form  of  the  displacements  and  rotations. 
Separating  the  strains  into  mid-plane  strains  and  curvatures 
gives, 


These  laminate  strains  are  strains  caused  by  the  deformation 
of  the  plate  only.  The  strains  due  to  the  hygrothermal 
stress  resultants  and  couples  can  be  included  by  using  the 
A,  B,  D matrices,  relating  these  quantities  to  strains.  As 
presented  in  ref.  [1]. 


It  should  be  noted  that  no  attempt  has  been  made  to 
include  the  transverse  normal  strains  in  the  stress  calcula 


With  the  laminate  strains  determined,  the  stresses 


within  each  lamina  are 


CHAPTER  FOUR 

NUMERICAL  ANALYSIS 

In  order  to  determine  the  importance  of  hygrothermal 
loads  in  analyzing  a laminated  composite  plate,  solutions  are 
presented  using  the  theory  developed  in  the  previous  two 
chapters.  The  examples  given  indicate  the  importance  of  hy- 
grothermal loads,  as  well  as  the  complexity  of  the  theory  re- 
quired to  obtain  satisfactory  engineering  solutions. 

To  efficiently  obtain  solutions  using  the  solution 
technique  given  in  Chapter  3,  a computer  program  was  de- 
veloped. The  program  is  listed  in  Appendix  D,  along  with  a 
brief  description  of  the  input  data  required  and  an  example 
output.  The  program  analyzes  rectangular  laminated  plates, 
symmetric  or  unsymmetric,  with  clamped,  simply- supported,  or 
clamped-simple  support  boundary  conditions.  Loads  included 
are  through  the  thickness  moisture  and  temperature  gradients 
and  uniform  lateral  pressures.  The  output  gives  the  deflec- 
tions, rotations,  stress  resultants,  stress  couples,  strains 
and  stresses  in  the  plate  coordinate  system.  Also  included 
are  options  to  perform  the  analysis  based  on  linear  plate 
theory  (LPT) , plate  theory  including  transverse  shear  de- 
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formation  (TSD) , and  plate  theory  including  transverse  shear 
deformation  plus  transverse  normal  strains  (TSD+TNS) . 


In  using  the  Rayleigh-Ritz  procedure  to  obtain  solu- 
tions to  the  governing  equations,  large  systems  of  linear 
equations  are  generated.  Only  the  digital  computer  can  effi- 
ciently solve  such  systems.  However,  an  important  part  of 
using  the  computer  as  an  analysis  tool  for  solving  equations 
is  verification  of  the  program  logic  or  code.  The  detailed 
results  of  the  many  example  cases  run  to  verify  the  program 
used  are  given  in  Appendix  A.  The  examples  analyzed  include 
classical  isotropic  plate  solutions  as  well  as  duplication  of 
many  of  the  results  published  by  Whitney  [5,  11,  12,  13]  for 
both  symmetric  and  unsymmetric  laminated  plates. 

During  program  verification,  a computational  error 
was  uncovered  in  the  single  surface  moisture  absorption  ex- 
ample presented  by  Pipes  et  al  in  reference  [1].  As  an  aid 
to  those  attempting  to  duplicate  these  particular  results, 
the  corrected  stress  plots  and  details  are  given  in  Appendix 
B. 

In  using  an  approximate  solution  procedure  with  a 
finite  number  of  terms  for  the  assumed  displacements  and 
rotations,  it  is  important  to  determine  the  convergence  to 
the  solution.  For  the  results  presented,  all  odd  terms  up 
through  m equal  to  seven  were  used.  The  convergence 
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of  the  solution  was  evaluated  by  calculating  the  percent 
change  in  the  lateral  defection  solution  t , contributed 

by  each  term  containing  the  highest  index. 

For  the  simply-supported  results  present,  using  terms 
through  m,n  = 7 , 

< 0. 1 Vo 

For  the  clamped  boundary  conditions, 

/<Vr  < 0.7  % 

The  following  material  property  data  from  reference 
[6]  for  the  T300/5208  graphite  epoxy  composite  were  used  in 
the  analysis. 


/ 

f 

* 

70°F 

350°F 

E1(psi)  \ 

21. OxlO6 

18.7xl06 

E2(psi) 

1.7xl06 

0.87xl06 

G12(psi) 

0.65xl06 

0.65xl06 

V12 

0.21 

0.21 

V21 

0.017 

0.010 

V23 

0.33 

0.33 

yT (in/in  PF) 

-0.21xl0~6 

-0.005x10 

YT (in/in  °F) 

16.0xl0“6 

21.8xl0“S 

(in/in  %WT) 
y 2 (in/in  %WT) 


ho(in) 


70°F 


6.67x10 


0.0055 


350°F 


6.67x10 


0.0055 


Temperature  dependent  material  properties  for  each 
lamina  were  based  on  linear  interpolation  of  the  above  data 
using  the  mean  lamina  temperature.  Linear  interpolation  is 
sufficient  for  most  problems  encountered  due  to  the  much 
higher  thermal  diffusivity  compared  to  moisture  diffusivity 
[14]. 


Comparison  of  Theories 


To  assess  the  importance  of  the  more  complex  plate 
theory  including  transverse  shear  deformation  and  transverse 
normal  strain,  the  lateral  deflections  of  several  T300/5208 
plates  subjected  only  to  a linear  moisture  gradient  were 
evaluated.  Figure  4.1  shows  the  center  deflection  of  four 
layer  symmetric  and  unsymmetric,  angle  ply,  simply- supported 
laminates  using  the  different  plate  theories  for  various 
A/B  ratios.  Figure  4.2  shows  the  deflections  of  the  same 
plate  for  different  A/H  ratios.  To  determine  if  thickness 
had  an  effect,  results  for  a 24  layer,  symmetric,  angle  ply 
plate  are  shown  in  Figure  4.3  for  A/B  ratios  and  Figure  4.4 
for  A/H  ratios. 
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to  Mq  * 1.0%  on  the  top  surface  at  Dt  = 0 . The  plate 
is  assumed  to  be  moisture  free  for  Dt  < 0 . Figures  4.7, 
4.8,  and  4.9  show  the  stresses  in  the  x,y  and  xy  direc- 
tions at  various  Dt's  for  a plate  without  in-plane  re- 
straints against  the  dilatational  expansions.  Figures  4.10, 
4.11,  and  4.12  show  the  stresses  with  in-plane  restraints. 

As  a basis  of  comparison.  Figure  4.13  shows  the  stress  dis- 
tributions through  the  plate  for  the  pressure  and  tempera- 
ture conditions  only. 


U 

u 


The  second  example  is  a simply- supported  4 layer, 

Q = + 45°  , angle  ply  laminate.  The  loading  is  a uniform 
pressure  of  1 psi,  a linear  temperature  gradient  of  from 
200°F  on  the  top  surface  to  100°F  on  the  bottom,  and  a 
single  surface  moisture  absorption  load  for  Mq  =1.0%  . 

The  plate  is  again  assumed  to  be  moisture  free  at  Dt  < 0 , 
and  unrestrained  in-plane.  Figures  4.14  and  4.15  show  the 
stress  distributions  for  a symmetric,  /451,-  45/2g,  laminate. 
Figures  4.16  and  4.17  show  the  stress  distributions  for  an 
unsymmetric,  /45^,-  45^,45^,-  45/  laminate.  Since  the 
stress  distribution  in  the  x and  y directions  are  ident- 
ical, they  are  plotted  together  in  one  figure.  For  compari- 
son, the  pressure  and  temperature  loads  only  stress  distribu- 
tions are  included  in  the  figures. 


Figures  4.18  through  4.21  show  the  stresses  for  a 

cross  ply  laminate  under  the  same  load 


4 layer,  9=0,  90° 


t 


B 
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conditions  as  the  previous  example.  In  this  case  the  shear 
stresses  are  zero  and  therefore  not  shown. 

Finally,  Figures  4.22  and  4.23  show  the  variation  in 
center  deflection  with  time  for  the  symmetric  angle  and  cross 
ply  four  layer  laminate  examples.  The  line  labeled  "P  + t 
Only"  represents  the  deflections  caused  by  the  pressure  and 
temperature  loads  only,  which  do  not  vary  with  time.  The 
straight  lines  drawn  between  data  points  do  not  represent 
the  actual  variation  in  deflections.  These  lines  are  in- 
cluded only  as  an  aid  to  visualize  the  varying  center  de- 
flection. 
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TABLE  4.2 

T300/5208  GRAPHITE  EPOXY  PLATE 
SIMPLY-SUPPORTED  ALL  EDGES 
4-LAYER  /45,  -45,  -45,  45/ 
A/B  =1  H = 0.022  IN 
M(z)  = 1.0  + 18.1818Z  (%WT) 


A/H 

(A=B) 

WCENTER  * 10 

(IN) 

LPT 

TSP 

TSD+TNS 

5.0 

0.00834 

0.00810 

0.00810 

10.0 

0.0335 

0.0329 

0.0329 

15.0 

0.0750 

0.0745 

0.0745 

25.0 

0.2084 

0.2079 

0.2079 

35.0 

0.4085 

0.4080 

0.4080 

50.0 

0.8337 

0.8332 

0.8331 

75.0 

1.876 

1.875 

1.875 

100.0 

3.335 

3.334 

3.334 

200.0 

13.34 

13.34 

13.34 

300.0 

30.01 

30.01 

30.01 

500.0 

83.37 

83.37 

83.37 

T300/5208  GRAPHITE  EPOXY  PL-irt 


SIMPLY-SUPPORTED  ALL  EDGES 


M(z)  « 1.0  + 18.1818Z  (*WT) 


"center  x l°3  <“> 

LPT 

TSD 

TSD+TNS 

10.75 

10.75 

10.74 

8.207 

8.207 

8.190 

6.386 

6.386 

6.367 

4.296 

4.295 

4.278 

3.431 

3.431 

3.418 

2.839 

2.839 

2.831 

2.397 

2.397 

2.392 

2.052 

2.052 

2.048 

1.188 

1.189 

1.187 

0.749 

0.749 

0.749 

. 
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TABLE  4.4 

T300/5208  GRAPHITE  EPOXY  PLATE 
SIMPLY-SUPPORTED  ALL  EDGES 
4-LAYER  /45/  -45,  45,  -45/ 
A/B  =1  H = 0.022  IN 
M(z)  = 1.0  + 18.1818Z  (%WT) 


WCENTER  10 


(IN) 


A/  n 
(A=B) 

LPT 

TSD 

TSD+TNS 

5.0 

0.0107 

0.0105 

0.0104 

10.0 

0.0430 

0.0425 

0.0423 

15.0 

0.0967 

0.0962 

0.0958 

25.0 

0.2685 

0.2680 

0.2669 

35.0 

0.5263 

0.5257 

0.5236 

50.0 

1.074 

1.073 

1.069 

75.0 

2.417 

2.416 

2.406 

100.0 

4.296 

4.295 

4.278 

200.0 

17.18 

17.18 

17.11 

300.0 

38.66 

38.66 

38.50 
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TABLE  4.5 

T300/5208  GRAPHITE  EPOXY  PLATE 
SIMPLY-SUPPORTED  ALL  EDGES 
24-LAYER  /45,  -45,  -45,  45/6s 
A/H  = 100  H = 0.132  IN 
M(z)  = 1.0  + 18.1818Z  (%WT) 


1 

A/B 

(A/H=100) 

"center  x 10  (IN» 

LPT 

TSD 

TSD+TNS 

0.33 

3.021 

o • (JzZ 

V • 0*4 

0.5 

2.298 

2.298 

2.298 

0.67 

1.784 

1.784 

1.784 

1.0 

1.201 

1.200 

1.200 

1.25 

0.9586 

0.9585 

0.9585 

1.5 

0.7932 

0.7932 

0.7932 

1.75 

0.6705 

0.6706 

0.6706 

2.0 

0.5744 

0.5746 

0.5746 

3.0 

0.3340 

0.3344 

0.3344 

4.0 

0.2110 

0.2113 

0.2113 

A/H 


— Figure  4.2 

l 
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T 300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY-SUPPORTED  ALL  EDGES 
24-LAYER  ANGLE  PLY 
A/H  =100  H = 0.132  IN 
R(Z)  * 1.0  + 18.1818aZ  (XWT) 

P = 0 PSI 


/43.— 43.-45.43/ 
LPT .TSD.TSD+TNS 


T300/-5208  GRAPHITE  EP0XY  PLATE 
SIMPLV-SUPP0RTED  ALL  EDGES 
24-LAYER  ANGLE  PLY 
A/B  = 1 H = 0.132  IN 
B(2)  - 1.0  + 18 .1818*2  (XUT) 

P » 0 PSI 


M/M 


Figure  4.5 


M0ISTURE  DISTRIBUTI0N 
SINGLE  SURFACE  ABS0RPTI0N 
6-LAYER  PLATE  H0»  0.0055  IN 


mm 
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M0ISTURE  DISTRIBUTI0N 
SINGLE  SURFACE  ABS0RPTI0N 
-LAYER  PLATE  H0»  0.0055  IN 


T300/5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
✓0*45#-45/*  UNRESTRAINED  IN-PLANE 
A » B A^H  - 100  Hq-  0.0059  IN 
Mn»  1.0  X P —2  PSI  T * 150  F 


DT-tXKT9 

OT-iXlO"4- 


OT«lXKr* 

OT-ixicr4 

DT-3XHT4 

DT-DdO"* 


T300/'5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
<'0.43.-45/',,  UNRESTRAINED  IN-PLANE 
A » B A^H  “ 100  0.0053  IN 

M0»  1.0  X P —2  PSI  T - 150  F 


DT-3X10"4 

0T-1X10"2 


DT-iXlO-3 

DT-1X10"4 


DT-1X1CT5 

OToixicr4 


DT-SXlfT4 

DT-1X10"* 


1300/5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
/0»45.-45/m  UNRESTRAINED  IN-PLANE 
A » B A/H  » 100  Ho  - O.OOSS  IN 
R0«  1.0  * P =-2  PSI  T = 150  F 


mmm 


T300/5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
✓0,43,-45/'*  RESTRAINED  IN-PLANE 
A « B A'H  ■ 100  0.0033  IN 

R0®  1.0  X P =*-2  PSI  T - 150  F 


DT  - 1X10* 


DT  * 1X10' 


T30Q/5208  GRAPHITE  EP0XV  PLATE 
CLAMPEO  ALL  EDGES 
/O, 43,-45/^  RESTRAINED  IN-PLANE 
A ® B A'H  • 100  fy®  0.0053  IN 
R,®  1.0  X P —2  PSI  T ® 150  F 


DT  - 1X10' 


IT  » 3X10' 


DT  - 1X10"* 


DT  - 1X10"3 


T300/5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
✓0# 45,-45/*  RESTRAINED  IN-PLANE 
A ® 0 A/H  «»  100  0.0035  IN 

R0-  1.0  X P =-2  PSI  T - 150  F 


0T  - 3X10' 
DT  - 1X10' 


DT  - 3X10 


DT  - 1X10' 


OT  - 1X10 
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T300/5208  GRAPHITE  EP0XY  PLATE 
CLAMPED  ALL  EDGES 
✓0,43,-434*  A - B A^H  « 100 
Hg*  0.0035  IN  P —2  PSI  T - 130  F 


T300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPP0RTED  ALL  EDGES 
'45 ,-45,-43 *45/  UNRESTRAINED  IN-PLANE 
A - B A/H  « 100  Ho*  0.0055  IN 
R«»  1.0*  P»-1PSI  T -150+4545 . 5*2  F 


1X10 


Figure  4.14 


DT  a 1X10 
OT  • 1X10' 


7300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLV  SUPP0RTED  ALL  EDGES 
/45 , -45 , -45 , 45/  UNRESTRAINED  IN-PLANE 
A = B A/H  = 100  Hq  = 0.0055  IN 
M0 - 1.0*  P=~1PSI  1=150+4545 .5*2  F 


Figure  4.15 


T300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPP0RTED  ALL  EDGES 
/-45 , 45 , -45 , 45/  UNRESTRAINED  IN-PLANE 
A = B A/H  » 100  0.0055  IN 

R0«  1.0*  P=-1PSI  T =150+4545 . 5*Z  F 


DT  - 1X10 
DT  «*  5X10 
0T  *»  1X10 
DT  a 1X10 


DT  » 1X10 
DT  « 5X10’ 
DT  » 1X10’ 


DT  » 1X10^ 
DT  « 1X10"" 
P + T ONLY 


DT  * 1X10"' 
DT  a 3X10" 
P + T ONLY 


DT  a 5X10 
DT  * 1X10 


T300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPPORTED  ALL  EDGES 
/-45# 45 #—45 #45/  UNRESTRAINED  IN-PLANE 
A = B A/H  = 100  Hq  = 0.0055  IN 
Ff0=  1.0*  P=-1PSI  T=1 50+4545 .5*2  F 


Figure  4.17 


T300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPP0RTED  ALL  EDGES 
✓0#90,90,0/'  UNRESTRAINED  IN-PLANE 
A = B A/H  * 100  0.0055  IN 

1.0*  P*-1PSI  T =150+4545 . 5*2  F 


P T ONLY, 


DT  * 1X10 
DT  - 3X10' 


DT  * 1X10 
P + T ONLY 


1X10 

5X10' 


DT  » 3X10 
DT  - 1X10' 
DT  - 1X10' 


T30G/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPPORTED  ALL  EDGES 
/O *90 *90*0/  UNRESTRAINED  IN-PLANE 
A » B A/H  « 100  0.0055  IN 

Rn»  1.0*  P*-1PSI  T *150+4545 . 5*2  F 


DT  » 1X10 
DT  » 1X10'1 


DT  - 1X10 


DT  - 1X10 


0 


T30Q/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPPORTED  ALL  EDGES 
'90,0.90*0/  UNRESTRAINED  IN-PLANE 
A = B A'H  » 100  fya  0.0055  IN 
R0a  1.0*  Pa-lPSI  T =150+4545 . 5*Z  F 


— jDT  - 1X10”®  rOT  » 5X10' 

P + T ONLY  /rDT  « lXlO^/rDT  - 1X10' 


P + T ONLY,  tDT  - 1X10 


T300/5203  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPP0RTED  ALL  EDGES 
V30»0»90#Q/  UNRESTRAINED  IN-PLANE 
A * B A/'H  * 100  H^ » 0.0055  IN 
PQ  « 1 .Or  P*-1PSI  T«150+454.5 . 5*2 


DT  - 1X10' 
DT  » 1X10' 
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T300/5208  GRAPHITE  EP0XY  PLATE 
SIMPLY  SUPPORTED  ALL  EDGES 
4-LAYER  SYM.  /45# -45 #-45 #45/ 

Re-  1.0*  P=-1PSI  T«150+4545.5*Z  F 
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CHAPTER  FIVE 

DISCUSSION  OF  RESULTS 

The  results  of  the  numerical  examples  outlined  in 
Chapter  Four  will  be  discussed.  The  several  points  of  in- 
terest will  be  described  in  the  following  two  sections  under 
the  importance  of  TNS  and  the  importance  of  hygrothermal 
loads  in  structural  analysis. 

Importance  of  TNS 

From  Figures  4.1  and  4.3  it  is  evident  the  effect  of 
both  TSD  and  TSD+TNS  on  the  center  deflection,  as  a function 
the  A/B  ratio,  is  small.  In  studying  these  figures,  it  is 
important  to  note  that  the  A/H  ratio  is  held  constant  at  100. 
Consequently,  only  B is  varied  to  obtain  various  A/B  ratios. 
If  the  center  deflection  were  non-dimensionalized  by  dividing 
by  B,  the  more  traditional  deflection  curve  with  a peak  at 
A/B  = 1,  would  be  obtained.  It  should  also  be  noted  that  the 
deflections  of  the  4-layer  unsymmetric  laminate  are  greater 
than  those  of  the  symmetric  laminate,  as  one  would  anticipate 
[11,  12,  13]. 

The  effect  on  deflection  of  the  width-to-thickness 
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ratio,  A/H,  is  also  small  when  comparing  the  three  theories 
for  A/H's  greater  than  100.  For  A/H  ratio's  less  the  100, 
Tables  4.2,  4.4  and  4.6  show  a slight  decrease  in  the  center 
deflection  when  including  TSD  compared  to  LPT. 

The  percent  decrease  in  the  deflection  using  TSD  as 
compared  to  LPT  for  the  three  example  plates  is  tabulated  be- 
low for  various  A/H  ratios. 


A/H 

% Decrease  in  WCE.NTER 

4-Layer 

Symmetric 

4-Layer 

Unsymmetric 

24-Layer 

Symmetric 

5 

3% 

3.8% 

3% 

25 

.2 

.2 

.3 

100 

0 

0 

0 

It  is  interesting  to  note  that  contrary  to  what  might 
be  expected  from  previous  efforts  with  uniform  lateral  loads 
[53 , the  magnitude  of  the  plate  deflections  decrease  when  in- 
cluding TSD  over  LPT.  The  reason  for  the  decrease  is  the 
type  of  loading.  For  hygrothermal  loadings,  the  dilatational 
expansions  are  satisfied  by  plate  curvatures  and  consequently 
deflections.  However  as  greatly  exaggerated  in  the  following 
sketches,  by  including  TSD,  part  of  the  expansions  are  accom- 
modated by  shear  deformations.  This  reduces  the  required 
curvature  change  and  thus  deflections,  necessary  to  satisfy 
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the  expansions. 


The  effect  of  the  TSD+TNS  theory  appears  to  be  neg- 
ligible. The  slight  decrease  in  the  plate  deflections  for 
the  unsymmetric  laminate  using  the  TSD+TNS  are  less  than  one 
percent,  and  not  considered  significant. 

In  comparing  the  results  for  the  4 and  24  layer  lam- 
inates, the  effect  of  laminate  thickness  is  also  negligible 
for  the  three  plate  theories. 


Importance  of  Hygrothermal  Loads 


The  stress  plots  shown  in  Figures  4.7  through  4.9  are 
for  the  6-layer,  /0,  45,  -45/2s  laminate  with  clamped  bound- 
ary conditions.  These  plots  illustrate  the  stresses  gen- 
erated by  pressure  acting  down  on  top  the  plate  along  with  a 
sudden  moisture  concentration  increase  on  the  top  surface. 


In  examining  Figures  4.7  and  4.8,  the  Sx  and  Sy 
on  the  top,  0=0°,  layer  are  compressive  as  expected. 
The  maximum  Sx  compressive  stress  of  approximately  9 ksi 
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occurs  in  the  fully  saturated  condition,  while  the  maximum 
Sy  stress  of  10  ksi  occurs  at  early  times.  Also  of  in- 
terest is  the  development  of  compressive  stresses  in  the 
bottom  0°  layer  in  the  saturated  condition. 

The  maximum  tensile  stress  of  approximately  5 ksi 
occurs  at  late  times  in  the  bottom,  0 = + 45°  , layer. 

The  shear  stress,  Sxy  , increases  with  increased 
moisture  concentration.  The  maximum  shear  stress  is  approxi- 
mately 12.5  ksi  with  large  shear  stress  gradients  between  the 
+ 45°  layers.  These  large  stress  gradients  imply  high  inter- 
laminar shear  stresses  as  noted  in  [1]. 

As  a comparison  of  the  hygrothermal  environment  ef- 
fects, Figure  4.13  shows  the  Sx  , Sy  and  Sxy  stress  dis- 
tributions for  the  same  plate  with  all  loads  except  moisture. 
Note  however  that  the  scale  for  these  curves  has  been  ex- 
panded when  making  comparisons.  The  most  obvious  differences 
from  Figure  4.13  are  the  increased  Sy  stress  in  the  0# 
layers  and  the  large  increase  in  Sxy  stresses  for  the  re- 
sults with  hygrothermal  effects  included  in  the  analysis. 

Figures  4.10  through  4.12  show  the  stress  distribu- 
tions for  the  identical  problem  as  before  but  with  in-plane 
expansions  restrained.  The  following  observations  can  be 
made  in  comparing  these  curves  to  those  of  the  previous  ex- 
ample . 
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1.  The  Sx  and  Sy  stresses  in  the  45°  layers  have 

shifted  negative  with  the  maximum  compressive  stresses  occur- 
ring at  the  fully  saturated  time. 

2.  The  Sx  and  Sy  stresses  have  not  changed  in 

the  0°  layers. 

3.  The  Sxy  shear  stresses  have  markedly  decreased 
in  the  middle  layers  with  a maximum  stress  of  only  5 ksi. 

The  general  shift  to  compressive  stresses  for  Sx 
and  Sy  in  the  middle  layers  is  expected  for  dilatational 
in-plane  restraints.  However,  the  maximum  compressive  stress 
is  only  slightly  higher  than  before  at  7 ksi. 

The  interesting  decrease  in  shear  stresses,  to  less 
than  half  those  for  the  unrestrained  case  shown  in  Figure 
4.9,  is  significant.  This  would  indicate  an  advantage  to 
designing  with  in-plane  restraints  to  lower  the  Sxy 
stresses  without  materially  increasing  Sx  and  Sy  . By 
performing  a buckling  analysis,  such  as  presented  by  Flaggs 
[6],  to  insure  plate  stability,  the  designer  could  reduce 
the  shear  stresses  and  shear  stress  gradients  by  providing 
in-plane  restraints.  This  would  also  reduce  the  possibility 
of  delamination  caused  by  the  interlaminar  shear  stresses. 

Table  5.1  illustrates  the  importance  of  the  hygro- 
thermal  loads  for  these  two  examples.  In  this  table  the  max- 
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TABLE  5.1 

COMPARISON  OF  MAXIMUM  STRESSES 
WITH  AND  WITHOUT  MOISTURE  LOADING 


COMPRESSION 

TENSION 

STRESS  RATIO 

shygro(ksi) 

LOCATION 

STRESS  RATIO 

shygro(ksi) 

LOCATION 

mm i 

z/Ho 

z/Ho 

Sx 

Fig. 4. 7 

4:1  - i-56 

3.0 

(0 "Layer) 

© 

00 

• 

O 

II 

-3.0 

(0 "Layer) 

Sy 

Fig. 4. 8 

^5  = 30‘6 

3.0 

(0 "Layer) 

13  * 3*47 

-2.0 

(45 "Layer) 

Sxy 

Fig. 4. 9 

-1-— l - 19.6 
• o 

-1.0 

(-45®Layer) 

33  = 9*6 

-2.0 

(45®Layer) 

Fig?4.10 

33  “ 1,36 

3.0 

(0 "Layer) 

rl-  i-0 

-3.0 

(0 "Layer) 

Sy 

Fig. 4. 11 

= 33.3 

3.0 

(0 "Layer) 

13  85  1*° 

-2.0 

(45®Layer) 

Sxy 

Fig. 4. 12 

^ - 7 5 
- .6  />5 

-1.0 

(-45°Layer) 

5.2  _ . n 
T3  ~ 4*° 

-2.0 

(45®Layer) 

0 

D 


I 

I 

I 

I 

I 
0 

II 


D 

D 


I 

I 


■ 
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imum  compressive  and  tensile  stresses  for  the  analysis  in- 
cluding hygro thermal  effects  are  compared  to  the  correspond- 
ing maximum  stress  for  the  same  laminate  without  hygrothermal 
loads . 

As  can  be  seen  in  Table  5.1  the  Sy  compressive 
stress  maximums  are  grossly  underpredicted  if  hygrothermal 
loads  are  not  included  in  the  analysis.  Another  area  in 
which  an  analysis  not  including  hygrothermal  effects  under- 
estimates the  stresses  is  in  the  middle  lamina  shear  stress- 
es, Sxy  . Only  for  the  Sx  maximum  tensile  stress  for  the 
unrestrained  laminate  does  the  analysis  not  including  moist- 
ure loading  predict  a higher  stress. 

Figures  4.14  through  4.17  show  the  stress  distribu- 
tions for  a simply-supported,  4-layer,  angle  ply  laminate 
subject  to  a lateral  pressure,  a linear  temperature  gradient, 
and  a moisture  gradient.  Figure  4.14  and  4.15  represent  the 
case  of  a symmetric  layup  pattern.  Figures  4.16  and  4.17 
are  distributions  for  the  unsymmetric  layup  pattern. 

The  Sx  and  Sy  stress  distributions  for  the  angle 
ply»  0 = + 45°  , laminates  are  identical  and  are  therefore 
plotted  together . Also  note  that  for  each  figure  the  curve 
labeled  "p+T  Only"  is  the  solution  for  the  identical  plate 
under  all  the  same  loadings  except  moisture. 
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Figure  4.14  for  the  symmetric  plate  shows  the  maximum 
Sx  and  Sy  stresses  occur  at  initial  moisture  introduction 
to  the  top  surface.  These  maximum  stresses  then  decrease  as 
moisture  is  absorbed  to  the  P+T  only  stress  values.  The 
shear  stress,  Sxy  , distribution  responds  in  the  opposite 
manner.  The  early  time  shear  stresses  are  approximately  the 
same  as  for  the  P+T  only  case.  However,  the  shear  stresses 
increase  and  reach  a maximum  as  the  plate  approaches  the 
fully  saturated  condition.  Again  note  as  with  the  previous 
example,  the  large  shear  stress  gradients  generated  at  the 
lamina  boundaries  with  moisture  absorption. 

The  stress  distributions  for  the  unsymmetric  laminate, 
starting  with  0 = - 45°  at  the  top  of  the  laminate,  are 
shown  in  Figures  4.16  and  4.17.  As  a general  observation, 
note  the  more  disjointed  pattern  of  stresses  at  the  lamina 
boundaries  for  the  unsymmetric  plate,  compared  to  the  stress 
distributions  for  P+T  only,  the  presence  of  the  moisture 
gradient  exaggerates  this  condition. 

The  Sx  and  Sy  stress  maximums  are  generally 
higher  for  all  layers  compared  to  the  symmetric  laminate.  In 
addition,  observe  that  the  Sx  and  Sy  stresses  do  not  ap- 
proach the  P+T  only  stresses  in  the  saturated  condition,  as 
for  the  symmetric  laminate. 


The  Sxy  shear  stress  distributions  in  Figure  4.17 
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respect  to  time.  In  these  figures  the  P+T  only  line  repre- 
sents the  deflection  caused  by  the  pressure  and  temperature 
gradient  only  and  does  not  vary  with  time.  As  would  be  ex- 
pected for  the  single  surface  moisture  absorption  the  de- 
flections decrease  initially  due  to  the  moment  resultants 
generated  by  the  moisture  gradient.  What  is  interesting  is 
that  as  the  plate  reaches  the  saturated  condition,  the  de- 
flections actual  increase  beyond  those  for  P+T  only.  The 
reason  for  this  reversal  is  the  moment  resultant  generated 
by  the  temperature  dependent  material  properties  of  the 
laminate.  These  properties  given  in  Chapter  Four  show  that 
the  stiffness  properties  of  the  composite  material  vary  with 
the  100°F  temperature  gradient.  By  looking  at  equation 
(2.26)  for  the  hygrothermal  moment  resultants,  one  sees  that 
even  though  the  moisture  distribution  is  constant,  the 
lamina  stiffnesses  are  not.  Thus,  a residual  moment  is  gen- 
erated which  adds  to  the  deflections  for  the  P+T  only  condi- 
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CHAPTER  SIX 

CONCLUSIONS 

An  analysis  of  generally  laminated  rectangular  com- 
posite material  plates  has  been  presented.  Loadings  in- 
cluded are  hygrothermal  gradients  and  lateral  pressures.  In 
developing  the  analysis  the  effects  of  transverse  shear  and 
transverse  normal  deformation  have  been  included.  The  re- 
sults presented  have  shown  the  importance  of  including  hygro- 
scopic loadings  in  designing  a composite  plate  structure,  as 
well  as,  the  degree  of  complexity  necessary  in  the  theory 
used  to  analyze  these  structures. 

The  results  of  the  parametric  study  of  the  importance 
of  transverse  shear  and  transverse  normal  deformation  have 
shown : 

1.  The  effect  of  transverse  normal  deformation  is 
negligible. 

2.  The  effect  of  transverse  shear  deformation  is 
less  pronounced  for  hygrothermal  loads  for  the  material 
studied  than  for  lateral  pressure  loads  presented  by  [5] . 


I 

I 

I 


3.  Transverse  shear  deformation  is  needed  to  analyze 
"thick"  (i.e.  A/H<100)  laminated  plates  which  are  subjected 
to  hygrothermal  and  pressure  loads. 

Stress  results  presented  show  that  hygrothermal  loads 
are  important  in  a detailed  analysis  of  a composite  plate. 

As  demonstrated  in  the  examples  given,  lamina  stresses  pre- 
dicted without  including  moisture  effects  may  be  in  error  by 
several  hundred  percent.  Based  on  the  diffusion  constants 
presented  in  [14]  and  the  plate  thickness,  the  maximum  lamina 
stresses  may  take  two  or  three  years  to  develop.  However, 
these  stresses  could  cause  premature  failure  if  not  accounted 
for  in  the  design  analysis. 

Hygrothermal  loads  have  also  been  shown  to  be  import- 
ant in  predicting  plate  deflections.  If  plate  deflections 
are  critical  to  a particular  design  hygrothermal  loads  must 
be  included  in  the  design  analysis. 
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APPENDIX  A 


VERIFICATION  OF  COMPUTER  PROGRAM 


The  computer  program  given  in  Appendix  D is  lengthy 
and  difficult  to  verify  by  step-by-step  comparisons  with 
hand  calculations.  Therefore,  to  verify  the  accuracy  of  the 
program  calculations,  several  example  problems  were  compared 
against  known  classical  isotropic  plate  theory  solutions  and 
solutions  published  in  the  literature  for  rectangular  lami- 
nated plates. 


The  following  solutions  for  a square  steel  plate 
under  various  load  conditions  were  evaluated  using  the  mate- 
rial properties, 

30.0  x 10^  psi 
0.3 

6.6  x 10-6  in/in  °F  , 

and  dimensions, 

A * B * 16  in 

H = 0.4  in  . 
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1 . Uniform  Pressure 

For  the  load  condition  of  a uniform  lateral  pressure 
of  1 psi,  the  center  deflection  is  compared  for  various 
boundary  conditions. 

Classical  TSD  + 

Theory  [36]  LPT  TSD  TNS 

Simply-Supported,  W^^io^  (in)  1.513  1.507  1.518  1.518 


Clamped,  WCTRxlo4(in) 


C lamped-Simple,  WCTRxlQ4(in)  7.157 


4.710  4.756  4.790  4.790 

7.154  7.197  7.197 


2..  Linear  Temperature  Gradient 

For  the  load  condition  of  a linear  temperature  gra- 
dient through  the  plate  thickness  of  T(z)  » 70  + lOOOz  °F 
with  Tppp  * 70°F  , the  center  deflections  for  various  bound* 
ary  conditions  af#  as  follows.  The  0late  dimensions  are 
A »• B * 100.0  in  and  h • 0.1  in  for  this  example. 

Classical 


Theory  (37] 

LPT 

TSD 

TSD+TNS 

MCTR(in)  6.19 

6.19 

6.18 

6.18 

0.0 

0.0 

0.0 

0.0 

Clamped,  wCTR(in) 


To  verify  the  stress  calculation,  the  following  prob- 
lems were  evaluated  using  the  same  material  properties  but 
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nated  plates  were  evaluated  using  LPT  and  TSD,  duplicating 
the  examples  presented  by  Whitney  [5].  The  following  mate- 
rial properties,  dimensions  and  load  data  were  used. 

4-layer  /0,  90,  90,  0/  A/B  - 1 


G12/E22  : 

= 0 

.6 

G23/E22 

= 0.5 

V12 

* 0.25 

Plx,y)  - 

Po 

sin  ££  sin 

A 

¥ ? 

B o 

■ 1 psi 

deflection  results 

are  non* 

-dimensionalized  using 

w' 

— 

p. 

X 

1C? 

Ell/E22 

s 

40 

Whitney 

W‘ 

111 

Sloan 

W* 

A/H 

= 

10 

LPT 

2.8 

2.8 

TSD 

4.4 

4.7 

A/H 

s 

40 

LPT 

2.8 

2.8 

TSD 

2.9 

2.9 

Ell/E22 

ss 

3 

A/H 

= 

5 

LPT 

18 

19 

TSD 

22 

25 

A/H 

25 

LPT 

18 

19 

TSD 

19 

19 

2.  Symmetric  and  Unsymmetric  Angle  Ply  Laminates 
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To  verify  the  program  for  evaluation  of  unsymmetric, 
as  well  as  symmetric  laminates,  several  examples,  both  sim- 
ply-supported and  clamped,  presented  by  Whitney  (12,  131 
were  duplicated.  The  load  used  in  the  exanqples  is  a uniform 
pressure,  PQ  = 1.0  psi  . 

The  deflection  results  for  the  symmetric  simply- sup- 
port plates  and  the  clamped  plates  were  non-dimensionalized 
using. 


AM  P. 


X to 


The  unsymmetric  deflection  results  were  also  non-dimens ion- 
alized  and  presented  as  a ratio  of  the  unsymmetric  compared 
to  the  symmetric  deflection. 


Simply-Supported : 


G12E22  * °‘5 


= 0.25  e = + 45« 


H = 0.011  in  A - B - 10.0  in 


Symmetric  Plate 


Whitney  112)  Sloan 


4-layer  /45,  -45,  -45,  45/  En/E22=40  0.235 

Ell/E22®10  0.875 


0.235 


0.875 
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Unsymmetric  Plate 

Whitney  [12] 

W* 

Sloan 

W' 

4- layer  /45,  -45,  45,  -45/ 

Su/E22-40 

1.2 

1.2 

eu/e22-io 

1.1 

1.1 

2r layer  /45,  -45/ 

Eu/E22-40 

3.1 

3.1 

Ell/E22-10 

2.0 

2.0 

Clamped: 

Elj/E22  * 40  G12//E22  * 0,5  v12  “ 0,25 

Q * + 45°  H = 0.011  in  A * B * 10.0  in 


Whitney  [ 13 ] 

Sloan 

W' 

W' 

Symmetric  Plate 

4-layer  /45,  -45,  -45,  45/ 

1.0 

1.0 

Unsymmetric  Plate 

4- layer  /45,  -45,  45,  -45/ 

1.2 

1.2 

2-layer  /45,  -45/ 

2.8 

2.7 

g I ; i ,t  • rt 

3.  Symmetric  and  Unsymmetric  Cross  Ply  Laminates 

In  the  same  matter,  several  cross  ply  symmetric  and 
unsymmetric  laminates  were  evaluated. 

' . 4 _ . . ; . ;■  . ;■  * ’ - .•  ....  * 

Simply-Supported : 

) ■ * n • •' • V.’  • 

Ell/E22  * 40  G!2/E22  = 1,0  v12  = °*25 


APPENDIX  B 


STRESS  PLOTS  FOR  SINGLE  SURFACE 
ABSORPTION  CASE  OF  REF.  [1] 


As  an  aid  to  those  attempting  to  duplicate  the  re- 
sults presented  in  [1] , the  correct  stress  plots  for  the 
single  surface  moisture  absorption  example  are  presented. 
The  error  described  was  made  in  the  numerical  calculations 
only.  The  equations  presented  in  [1]  are  believed  to  be 
correct. 


The  stresses  within  each  lamina  can  be  calculated 
according  to  equation  (27)  of  [1],  or 

** “ftw-)]  (B.D 

‘.j- '.*,<• 

The  stress  computations  of  [1]  omitted  the  z.X*(>0  term. 
This  omission  only  effected  the  single  surface  absorption 
results  because  of  the  nonzero  curvatures.  The  other  ex- 
amples presented  had  symmetric  loading  and  therefore  zero 
curvatures,  Xj 

Another  difference  between  the  presented  plots  and 


those  of  [1]  is  in  the  sign  of  the  shear  stresses.  The  fol- 
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T300/5208  GRAPHITE  EP0XY  PLATE 
RECTANGULAR  PLATE  ELEMENT 
/O# 45 .*-45/' a*  = 0.0055  IN  = 1.4% 


T300/5208  GRAPHITE  EP0XY  PLATE 
RECTANGULAR  PLATE  ELEMENT 
/0#45*-45/«»  H,  * 0.0055  IN  = 1.4% 


IT  * 1X10 


7 ■ 1X10' 


DT  - iXiO' 


•oo  laoo 


(KSI ) 


DT  - 1X10^ 
DT  - 1X10“* 


DT  - 1X10*"4 


DT  - 1X10“* 


T300/5208  GRAPHITE  EP0XY  PLATE 
RECTANGULAR  PLATE  ELEMENT 
✓G, 45,-45/2,  Ho  = 0.0055  IN  1.4* 


Figure  B. 3 
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lowing  stress  plots  are  based  on  the  definition  of  6 as 
shown  in  Figure  2.1.  Using  this  definition  the  sign  of  the 
shear  stresses  is  opposite  to  that  presented  in  [1]. 

Figures  B.l,  B.2  and  B.3  show  the  correct  stress 
distributions  for  the  single  surface  absorption  example. 
Figures  B.l  and  B.2  of  the  x and  y stress  distributions, 
respectively,  show  a much  larger  stress  gradient  particularly 
in  the  middle  layers.  The  shear  stress  distributions  of 
Figure  B.3  also  indicate  a more  bending  type  stress  gradient 
through  the  center  layers.  It  is  important  to  note  the  shear 
stress  in  the  6*0°  layers  is  not  zero  for  early  times . 
This  nonzero  shear  is  generated  by  the  unsymmetric  moisture 
distribution  in  the  +45®  layers  even  though  the  layup  pat- 
tern is  symmetric.  For  the  fully  saturated  condition, 

_2 

Dt  = 1 x io  , the  presented  curves  match  those  of  [1]  ex- 
cept for  the  sign  of  the  shear  stresses. 
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APPENDIX  C 

EVALUATION  OF  THE 

CHARACTERISTIC  BEAM  FUNCTION  INTEGRALS 

The  presented  evaluation  of  the  general  integral 
forms  of  the  characteristic  beam  functions  are  provided 
to  assist  those  attempting  to  duplicate  this  work.  The 
characteristic  beam  functions  used  in  Chapter  Three  are 
expressed  in  the  following  shorthand  notation. 


1:  <fc  w = (f-  - i)  + $ cosK /“*(■§•  - 

i (t  - 4)  - c“*k  A (£  - i) 

9;(Z)  » (§  . -V  7*  /**(%-  ± ) 

j 

e>)  = 

cc . (z)  = cos  ^LA_S 

^ C 

SSi(x)  = sii^  Jjk* 


U 

[] 

0 


where  z = dummy  coordinate  direction 

c = length  of  plate  along  z direction 
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dummy  mode  index*  odd  only. 


As  given  in  equations  3.11  and  3.12* 


or 


_ Si  A (fi/Z  ) 


and 


A 


are  the  solutions  to 


+o.v\  (/fc/x)  i*  +<xvnV\ 


(C  .3) 


(C.4) 


^ ' ' A*3' 


(C.5) 


The  following  forms  of  expressions  occur  often  in 
the  evaluations  and  are  abbreviated  as 
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0 


01: 
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+ smOW.jVfc 


(Pi+Pri/z 


NN-.  = ft «®  & sUv>^.  + ftj  si'^tjL 
Pf-j  - coS.Uf£  _ fiCo »jfdi  si»k ^ 


(C.10) 


(C.ll) 


(C.12) 


(C.13) 


(C.14) 


(C.1S) 


(C.16) 


(C.17) 


(C.18) 


(C.19) 


beam  functions  across  the  domain  of  interest  are  needed  in 


applying  the  Rayleigh-Ritz  solution  technique.  Thus,  the 
following  general  integral  forms  are  useful  in  evaluating 
the  necessary  terms 
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- cHH«  +?« 


1-1*  c LL.ri  > 

J^rW*  = 2 - ^ CCr  . 


v>s  (C.25) 


'TteS  (C. 26) 


n 


C LL 


rs  » 


ir*s  (C. 27) 


fssr<f>tJz  - c Si* as  E E + PF  , ^f(c.28) 

© **  /C  * ** 

A-4i«J*.  = CI,;'9t  E6„-  g^*'"?  pp  ~?<C.29) 
* ft,*  + W)‘  ^ 


)“■& 4l  =-CS«^  2fl<LS(~ 


/l (V  ^*(C.30) 
,N(r^  ^ * **8 


{We;  i.-cwi  G-G  , 

/*S  * V^W) 


APPENDIX  D 


RECTANGULAR  COMPOSITE  MATERIAL  PLATE 
COMPUTER  PROGRAM 


A computer  program  for  the  analysis  of  laminated 
composite  material  plates  has  been  written.  The  program, 
LAMINATE,  performs  the  numerical  computations  for  analyzing 
rectangular  plates  using  the  solution  technique  given  in 
Chapter  3.  LAMINATE  will  handle  plates  with  all  edges  sim- 
ply-supported, or  all  edges  clamped,  or  with  two  opposite 
edges  clamped  and  the  other  two  edges  simply-supported.  The 
program  will  analyze  plates  subject  to  lateral  pressures  and 
with  temperature  and  moisture  distributions  through  the 
thickness.  Composite  material  plates  may  be  up  to  32  layers 
thick  of  either  symmetric  or  unsymmetric  layup  patterns. 


Description  of  Input  Data 

The  input  data  required  to  run  problems  using  the 
LAMINATE  computer  code  is  assembled  as  a data  card  deck.  A 
data  card  deck  consists  of  a card  type  A and  any  number  of 
data  card  sets  (card  types  B through  J) . One  data  card  set 
is  included  for  each  case  of  a run. 
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The  first  card  of  a data  card  deck  must  be  a type  A 
card,  as  this  card  inputs  the  number  of  data  card  sets  to  be 
read.  After  card  type  A,  NCASE,  data  card  sets  (card  types 
B through  J)  must  be  included. 

Mot  every  data  card  set  needs  to  contain  all  the  card 
types  B through  J.  Depending  on  the  control  options  select- 
ed, some  of  the  card  types  may  be  omitted.  These  control  op- 
tions allow  data  which  was  inputted  in  the  previous  case  to 
be  omitted  from  subsequent  data  cards  sets,  if  desired. 

All  cards  of  card  types  B,  C,  E,  F,  G,  and  H must  be 
included  in  each  data  card  set.  Cards  of  card  types  D and  I 
may  be  omitted  if  it  is  desired  to  use  values  from  the  pre- 
vious case. 

In  the  following  section  the  detailed  input  instruc- 
tions for  each  card  type  are  given.  If  there  is  more  than 
one  card  for  a particular  card  type#  the  cards  are  numbered 
in  the  order  in  which  they  should  appear  in  the  data  card 
set  (i.e.  Dl,  D2,  D3,  etc.).  The  information  given  for  each 
card  of  each  card  type  consists  of  the  variable  name,  the  in- 
put format,  and  a description  of  the  input  and  units  of  mea- 
sure required. 

All  input  for  an  integer  (I)  format  field  must  be 
right  justified.  All  input  for  exponential  (E)  format  must 
have  the  exponent  right  justified  in  the  field. 
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INPUT  DATA 

Card  Type  A (one  card)  must  be  included  only  once  £or  each 
computer  run  and  be  the  first  card  of  the  data  card  deck. 

Format  Description 

15  Number  of  B through  J data  card  sets 

for  this  run. 

Card  Type  B (one  card)  must  be  included  in  each  data  card 

set. 

Variable  Format  Description 

COM ( i ) * 18A4  Comment  up  to  72  characters  long;  any 

alphanumeric  character  except/  (;), 
semicolon . 

Card  Type  C (two  cards)  must  be  included  each  data  card  set. 

Description 

1 Simply- supported  all  edges 

2 Clamped  all  edges 

3 Clamped  x-edges,  simply- supported 
y-edges 

4 Rectangular  plate  element  per 
(I],  no  boundary  conditions 

indicates  subscripted  variable. 


Card  Cl 

Variable  Format 

IFUNCT  15 


Variable 

NCASE 


Variable  Format 


KTS  I5 


i 


i 


'KTNS  X5 

I 


KTBMP  15 


i 
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Description 

5 Calculate  and  print  the  A,  B,  o 
matrices  only 

- 1 Transverse  shear  included  using 
the  distribution  f(z)  * 

U-(f)2J 

2 Transverse  shear  included  using 
the  distribution  f(z)  = 

f 

3 Transverse  shear  including -with- 
out a weighting  distribution 

4 Simulate  the  neglect  of  trans- 
verse shear  deformation 

- 0 Neglect  transverse  normal  strains 
• 1 include  transverse  normal  strains 

- 0 Ho  temperature  dependent  material 

properties 

1 Linearly  interpolate  material 
properties  based  on  the  mean  lam- 
ina temperature  (KMATP  must  equal 
2) 

Number  of  points  at  which  a detailed 
deflection  and/or  stress  printout 
is  desired  (number  of  Card  Type 
J s included  in  data  card  aet) 


109 


K 
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Variable  Format  Description 

KHYGRO  15  = 0 Neglect  all  moisture  and  tempera- 

ture loading 

* 1 Polynomial  temperature  distribu- 
tion  only 

■ 2 Polynomial  moisture  distribution 
only 

= 3 Both  polynomial  temperature  and 
moisture  distributions 

KPRESS  15  =1  Uniform  lateral  pressure 

= 2 Whitney  pressure  distribution 

P(x,y)  = PQ  sin  sin  ^ 

= 3 No  lateral  pressure  loadings 
KMATP  15  I KMATP | is  the  number  of  material 

property  data  sets,  card  type  D's. 
If  KMATP  is  positive,  KMATP  is 
the  number  of  Card  Type  D sets 
inputted. 

If  KMATP  is  negative,  | KMATP | is  the 
number  of  Card  Type  D's  used  in 
the  previous  case,  or  | KMATP | is 
the  number  of  laminate  layers  if 
in  the  previous  case  material 
properties  were  linearly  inter- 
polated (KTEMP-1) . 


Variable  Format 

ITERMN  15 


KMOIST  15 


KSYM  15 

KCROSS  15 

Card  C2 

Variable  Format 

15 

15 


KPRINT 


KPQ 


Description 

Maximum  odd  index  used  in  the  series 
expansion  of  the  assumed  dis- 
placement and  rotation  functions 
(ITERMN<9) 

= 0 Moisture  distribution  not  per 
diffusion  equation  solution 
= 1 Moisture  distribution  per  diffu- 
sion equation  solution  of  [1] . 
Note  KMOIST  = 1 will  override 
polynomial  moisture  distribution 
= 0 Unsymmetric  laminate 
= 1 Symmetric  laminate  (sets  = 

0) 

= 0 Angle  ply  or  mixed  laminate 
= 1 cross  ply  laminate 

Description 

= 0 Print  material  property  data 

■ 1 Suppress  all  material  property 

data  printout 

■ 0 Print  Q matrix  for  all  layers 

(see  Card  Type  I to  suppress 


■ 


nry 


Li 

u 

d 

Uj 

ii 

ii 

o 

li 

o 


printout  for  specific  layers.) 

1 Suppress  Q printout  for  all  layers 


Ill 


Variable 


Format 


Description 


KPQB 


KCSPR 


= 0 Print  Q matrix  for  all  layers 
(see  Card  Type  1 to  suppress 
printout  for  specific  layers) 

= 1 Suppress  Q printout  for  all 
layers 

= 0 Suppress  printout  of  simultaneous 
equations  coefficient  and  source 
matrices 

= 1 Print  simultaneous  equations  co- 
efficient and  source  matrix. 


KINTP 


(Note  number  of  lines  printed  is 
large  for  ITERMN>5.) 

= 0 Suppress  integration  subroutine 


check  run 


1 Perform  integration  subroutine 
check 


Card  Type  D (four  cards) . The  Card  Type  D cards  are  only  re- 
quired if  KMATP  >1.  If  KMATP  (card  Cl)  is  greater  than  or 
equal  to  one,  the  KMATP  is  the  number  of  Card  Type  D sets 
(i.e.  Dlr  D2 r D3 , D4)  inputted  for  in  case. 


Card  D1 


Variable 


Format 


Description 


TMAT(i) 


F10.3  Temperature  (°F)  for  the  material 

properties  in  this  Card  Type  D set 
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Variable  Format  Description 

EEI(i,j)  3E10.4  Elastic  moduli  in  E^,  E22'  E33 

order  (psi) 

Card  D2 

Variable  Format  Description 

GGl(i,j)  3E10.4  Shear  moduli  in  G12'  G13'  G23  order 

(psi) 


Card  D3 

Variable  Format  Description 

XNU(i,j)  6F10.3  Poisson  ratios  for  material  in-order 

V12'  v21'  v13 ' v31'  v23'  V32 

Card  D4 

Variable  Format  Description 

ALPHA (i,j)  3E10.4  Coefficients  of  thermal  expansion  in 

the  order  y1x»  Y22*  Y33  (in/in 

•F) 

BETA ( i / j ) 3E10.4  Coefficients  of  hygroscopic  expan- 

sion in  the  order  $11#  e22,  333 
(in/in  %WT) 

Card  Type  E (one  card)  must  be  included  in  each  data  card 
set. 


u 

Li 

u 

D 

D 

D 

D 

D 

0 
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li 

11 

0 

11 

0 

D 

11 

0 

u 

D 

D 

0 

D 

0 

e 
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KLAYER  is  negative , the  laminate 
of  | KLAYER | layers  used  in  the 
previous  example  will  be  used. 
P10.4  Plate  length  in  X-direction  (in) 

F10.4  Plate  length  in  Y-direction  (in) 

F10.4  Total  thickness  of  plate  (in) 

(one  card)  must  be  included  in  each  data  card 


Variable 

Format 

Description 

PRESS 

F10.4 

Lateral  pressure  (psi) , positive  act- 
ing down  on  plate 

TREF 

F10. 4 

Reference  temperature  of  plate  prior 

to  loading  (°F) 

XMREF 

F10. 4 

Reference  moisture  concentration  of 

plate  prior  to  loading  (%WT) 

XMS 

F10.4 

Moisture  concentration  for  diffusion 

solution  moisture  distribution 

(%WT) . Only  included  if 

KMOIST  > 0 on  Card  Cl. 

DT 

F10.4 

Diffusion  time  parameter  for  diffu- 

sion  equation  moisture  distribu- 
tions (in2) . Only  included  if 


AX 

BY 

H TOTAL 


set. 


KMOIST  > 0 on  Card  Cl 
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set. 


Variable 


TD(i) 


set. 


Variable 


HM(i) 


(one  card)  must  be  included  in  each  data  card 


Format 


Description 


5E12.3  Coefficients  for  the  temperature  dis- 
tribution through  the  thickness 
of  the  plate  in  the  form, 

Tex)  » Tt>(0  * Tb(iWZ  tTb(j)#Zx 
t Tt>(*) # za  * Tb^  a z* 

(one  card)  must  be  included  in  each  data  card 


Format 


5E12.3 


Description 

Coefficients  for  the  moisture  con- 
centration distribution  through 
the  plate  thickness  in  the  form, 

s «•  Wr\(i>Z  ♦ Mfn(j)az*" 


Card  Type  I (|KLAYER|  cards).  Card  Type  I's  are  only  in- 
cluded if  KLAYER  on  Card  Type  E is  greater  than  zero.  If 
KLAYER  > , then  KLAYER  cards  will  be  inputted,  one  for 

each  in  the  order  of  the  stacking  in  the  laminate 

at  the  bottom  (z— HTOTAL/2)  of  the  laminate. 


ciable 
THETA (i) 


F10.4 


Description 

Angle  from  the  X to  1 axis  for  this 
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Variable 


THICK (i) 


KMAT(i) 


Format 


F10.4 


Description 


lamina,  in  degrees 

Thickness  of  this  lamina  of  the  lam- 
inate (in) 

Material  property  (Cards  Dl  through 
D4)  data  set  for  this  lamina. 


KPROP(i) 


If  KTEMP  > 0, 


KMAT  must 


equal  the  layer  number  starting 
with  1 for  the  bottom  layer. 

0 Print  material  properties  for 


this  lamina 


KPQBL(i) 


= 1 Suppress  material  properties 
printout  for  this  lamina 
= 0 Print  stiffness  properties  for 
this  lamina 

* 1 Suppress  print  of  stiffness  prop- 
erties for  this  lamina 


Card  Type  J (KDET  cards).  KDET  Card  Type  J's  are  inputted 


in  each  data  card  set. 


Variable 


Format 


Description 


F10.3 


F10.3 


X-position  for  evaluation  of  the  de- 
flections and  stresses  (in) 
Y-position  for  evaluation  of  the  de- 
flections and  stresses  (in) 


Variable 


KSTRES 


Format 


KPRINV 


KPREL 


Description 

* 0 No  stress  calculation  for  this 

position  X,Y 

> 0 KSTRES  is  the  number  of  equally 
spaced  points  through  the  thick- 
ness for  each  lamina 
= 0 Suppress  printing  of  the  inverted 
A,  B#  D matrix 

* 1 Print  inverted  A,  B,  D matrix 

* 0 Suppress  printing  the  plate 

strains  (e)  and  curvatures  (k) 

= 1 Print  plate  strains  (e)  and 
curvatures  (tc) 

= 0 Suppress  printing  lamina  strains 
= 1 Print  lamina  strains  at  each 

stress  evaluation  location  through 
the  lamina 

= 0 No  in-plane  dilatational  plate 
restraints 

* 1 Restrain  in-plane  dilatational 

strains 


Repeat  Card  Types  B through  J as  required,  NCASE-1  times. 


KPREK 


KIPE 


>le  Probl< 


The  input  data  required  and  output  printout  for  the 


0 

0 

0 

D 
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following  example  problem  is  included  to  illustrate  the  use 
of  the  LAMINATE  program. 

The  example  problem  will  be  a 4- layer,  angle  ply# 
T300/5208  Graphite  Epoxy  laminated  plate.  All  edges  of  the 
plate  will  be  simply- supported.  The  plate  layup  pattern  is 
/45,  -45/s.  The  length  in  the  x and  y directions  is 
10.0  inches.  The  total  plate  thickness  is  0.022  inches. 
In-plane  dilatational  strains  will  be  unrestrained.  The 
plate  loading  will  be  a uniform  pressure#  PQ  * 0.1  psi  # a 
linear  temperature  gradient#  T(z)  = 150  + 4545.5  * z (°F) 
and  a moisture  distribution  based  on  the  single  surface  dif- 
fusion equation  solution  with  MQ  **  1.0%  WT  and  Dt  = 

1 x 10“5in2. 

The  input  data  card  data  for  the  example  is  shown  in 


Figure  D.l.  The  printout  (10  pages)  of  the  example  is  in- 
cluded after  Figure  0.1. 


